Structures and Stabilities of CaO and MgO Clusters and Cluster Ions: An alternative 

interpretation of the experimental mass spectra. 



O 

o 
o 

(N 
(N 



C/3 

s 

o 

• I— I 

^ ■ 

> 

^ ■ 

^ : 

o . 

o 

o , 
o , 

o 

^ : 



X 



Andres Aguado* and Jose M. Lopez 
Departamento de Fisica Teorica, Universidad de Valladolid, Valladolid 47011, Spain 

The structures and relative stabilities of doubly-charged nonstoichiometric (CaO)nCa^+ (n=l- 
29) cluster ions and of neutral stoichiometric (MgO)n and (CaO)n (n=3,6,9,12,15,18) clusters are 
studied through ab initio Perturbed Ion plus polarization calculations. The large coordination- 
dependent polarizabilities of oxide anions favor the formation of surface sites, making the critical 
cluster size where anions with bulk coordination first appear larger than that found in the related 
case of alkali halides. Thus, we show that there are substantial structural differences between alkali 
halide and alkaline-earth oxide cluster ions, contrary to what is suggested by the similarities in 
the experimental mass spectra. An alternative interpretation of the magic numbers for the case of 
oxides is proposed, which involves an explicit consideration of isomer structures different from the 
ground states. A comparison with the previously studied (MgO)nMg'^^ cluster ions shows that the 
emergence of bulklike structural properties with size is slower for calcium oxide. Nevertheless, the 
structures of the doubly charged clusters are rather similar for the two materials. On the contrary, 
the study of the neutrals reveals interesting structural differences between MgO and CaO, similar 
to those found in the case of alkali halides. 

PACS numbers: 36.40.C; 61.46. -fw; 61.50.Lt; 61.60. -fm; 79.60.Eq 



I. INTRODUCTION 

Small clusters are of great interest both to the phys- 
ical and chemical communities because of their numer- 
ous potential applications (for example, in nanoelectron- 
ics or catalysis), and also because one can gain impor- 
tant insight into the evolution from atomic and molec- 
ular properties to bulk and surface properties. To have 
a knowledge of the structures adopted by the clusters 
is of paramount importance, as many interesting cluster 
properties are largely determined by them. From the the- 
oretical side, finding the lowest energy structure for each 
cluster is a complicated matter, because the number of 
isomers increases exponentially with cluster size. Other 
reasons are that one has to treat bulklike and surfacelike 
ions on an equal footing, and that the number of ions 
to be explicitely considered in a cluster is larger than 
in a bulk or surface study, where symmetry restrictions 
impose a number of useful atomic equivalences. From 
the experimental side, the problem is so difHcult that 
during the approximately 30 years of intensive cluster 
research, the main source of structural information has 
been theory. Very recently, experimental techniques like 
electron diffraction from txajBped clustersEI or measure- 
ments of cluster mobilitiegJlS have been succesfully ap- 
plied to study the structures of covalent and ionic clus- 
ters. Photoelectron spectroscopy has also been applied 
to study- isomerization transitions in small alkali halide 
clusters,!] and measurements of ionization potential ta 
detect structural transitions in barium oxide clusters.^ 
At the moment, however, these techniques need parallel 
theoretical calculations to make a definite assignment of 
the observed diffraction pattern, mobility or ionization 
potential to a specific isomer geometry. 



A large amount of theoretical work has been devoted 
to metallic, semiconductor and noble gas clusters. The 
work on ionic materials has been centered mostly in the 
family of alkali halides, and studies of metal oxide clus- 
ters have been comparatively scarce, despite their impor- 
tance in many branches of surface phy^|iGS, like heteroge- 
neous catalysis or corrotion. SaundersuiiJ reported mass 
spectra and collision induced fragmentation data for stoi- 
chiometric4MgC))n ^^'^ (CaO)^ cluster ions, Martin and 
BergmannlliJ published mass spectra Cj£-(GaO)„Ca^+ clus- 
ter ions, and Ziemann and Castlemar£3~1l3 performed ex- 
perimental measurements of several singly- and doubly- 
ionized cluster ions of MgO and CaO by using laser- 
ionization time-of-flight mass spectrometry. Theoretical 
calculations have been performed at different levels of ac- 
curacy: simple ionic models based on phenomenologinal 
pair potentials were used by Ziemann and CastlemartJ^tj 
to explain the_.global trends found in their experi- 
ments; WilsonEj has studied neutral (MgOLj (n<30) 
clusters by using a compressible-ion modeltZI that_itti 
eludes coordination-dependent oxide polarizabilities jlS'ta 
semiempirical tight-binding calculations fo r Me O clusters 
were reported by Moukouri and NogueraJ2!jl22l finally, ab 
initio calculations on stoichiometric MgO-elusters have 
been presented, recently|— by Recio et aL,E3'E3 Malliavin 
and Coudray,Eil Li et a/,t3 and de la Puente et aljEB and 
calculations on stoichiometric (Li2ClLi„ clusters have been 
reported by Finocchi and Noguera.EZl Regacding the non- 
stoichiometric cluster ions, Aguado et alx^ have studied 
the structures and stabilities of (MgO)„Mg^+. 

Trying to find an interpretation oLthe obtained mass 
spectra, Ziemann and Castlemaro performed some 
simple pair potential calculations of the structures of 
(MgO)„Mg^+ cluster ions by using a rigid ion model. 
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The conclusion of those calculations was that the magic 
numbers can be explained in terms of highly compact 
structures that can only be obtained for certain cluster 
sizes, an interpretation very similar j-tn -that found in the 
closelw-jrelated case of alkali halides-cHtll In our previous 
work,ca we showed that the structures of (MgO)„Mg^+ 
cluster ions were quite different from those of alkali 
halides. Specifically, the influence of the large and coor- 
dination dependent polarizabilities of oxide anions (not 
included in the rigid ion model) favors the formation of 
surface oxide sites, and thus structures with bulk ox- 
ide anions (coordination 6) are not energetically com- 
petitive until large values of the number of molecules 
n are attained. For example, a highly compact 3x3x3 
cube (where the notation denotes the number of atoms 
along three perpendicular edges) is not the ground state 
of (MgO)i3Mg^+. Nevertheless, the agreement between 
the magic numbers obtained through an examination of 
the stabilities of the clusters against the loss of an MgO 
molecule and the experimental ones is complete. It is just 
the interpretation of them in terms of structures that is 
different, that is, model-dependent. It is interesting to 
study a similar system like calcium oxide in order to as- 
sess whether those trends are a general feature of alkaline- 
earth oxidCpClusters or not. Moreover, the experiments 
of SaunderELIilj suggest interesting structural differences 
between both materials, as the main fragments observed 
after collisions with inert gas ions were (MgO)3 in one 
case and (CaO)2 in thepatiier, and the mass spectra of 
Ziemann and CastlemanU'EJ show different stabilities in 
the small size regime (magic numbers at n=5,8,ll for 
(MgO)„Mg2+ and at n=5,7,9,ll for (CaO)„Ca2+), pro- 
viding further motivation for our study. From the theo- 
retical point of view, Ca^+ is larger than Mg^+ , so we can 
expect ionic size packing effects to play an important role 
in determining structural differences. Besides, Ca^"*" has 
a polarizability approximately 6 times larger than Mg^+ , 
and the polarizabilities of the oxide anions are also larger 
in CaO because the bonding is weaker than in MgO. 

In this work we present the results of an extensive and 
systematic study of (CaO)„Ca^+ cluster ions with n up 
to 29, and of (MgO)„ and (CaO)„ neutral stoichiometric 
clusters with n=3,6,9,12,15,18. The rest of the paper is 
organized as follows: in Section II we give a brief resume 
of the theoretical model employed, as full an jfiiposition 
has been reported already in previous works.Ej The re- 
sults are presented in Section III, and the main conclu- 
sions to be extracted from our study in Section IV. 

II. THE AIPI MODEL AND POLARIZATION 
CORRECTIONS 

The theoreticalj— foundation of the ab initio per- 
turbed ion -modelo lies in the theory of electronic 
separability,L3'Efl and its practical implementation in the 
Hartree-Fock-|^HF) version of the theory of electronic 
separabilityEaO Very briefly, the HF equations of the 



cluster are solved stepwise, by breaking the cluster wave 
function into local group functions (ionic in nature in 
our case). In each iteration, the total energy is min- 
imized with respect to variations of the electron den- 
sity localized in a given ion, with the electron densities 
of the other ions kept frozen. In the subsequent itera- 
tions each frozen ion assumes the role of japnfrozen ion. 
When the self-consistent process finishes,E^ the outputs 
are the total cluster energy and a set of localized wave 
functions, one for each geometrically nonequivalent ion of 
the cluster. These localized cluster-consistent ionic wave 
functions are then used to estimate the intraatomic cor- 
relation energy correc.tinn. through dementi's Coulomb- 
Hartree-Fock method.tllS The large multi-zeta basis sets 
of Clementi and RoettiEa are used for the description of 
the ions. At this respect, our optimizations have been 
performed using basis sets (5s4p) for Mg^+ and (5s5p) 
for , respectively. Inclusion of diffuse basis functions 
has been checked and shown unnecessary. One important 
advantage coming from the localized nature of the model 
is the linear scaling of the computational effort with the 
number of atoms in the cluster. This has allowed us to 
study clusters with as many as 59 atoms at a reasonable 
computational cost. ■— ■ 
In our previous work on alkaline-earth oxide clusters,cS 
we concluded that the aiPI model is equivalent to a first- 
princioles version of the semiempirical breathing shell 
model.E3 The binding energy of the cluster can be written 
as a sum of deformation and interaction terms 

Ebind = ^ E^ind = ^i^def + 2^mt)- (1) 
R R 

where the sum runs over all ions in the cluster. The 
interaction energy term is of the form 

T^R _ pRS _ / pRS I j^RS I j^RS , j^RS \ 

^int — 2^ ^int ~ 2^ \^ class + ^nc + + ^ overlap )^ 

S^R Sj^R 

(2) 

where the different energy contributions are: the clas- 
sical electrostatic interaction energy between point-like 
ions; the correction to this energy due to the finite exten- 
sion of the ionic wave functions; the exchange interaction 
energy between the electrons of ion R and those of the 
other ions in |tlje cluster; and the overlap repulsive energy 
contributioncJ The deformation energy term E^^^ is the 
self-energy of the ion R. It is an intrinsically quantum- 
mechanical many-body term that accounts for the energy 
change associated to the compression of the ionic wave 
functions upon cluster formation, and incorporates the 
correlation contribution to the binding energy. As the 
model assumes, for computational simplicity, that the 
ion densities have spherical symmetry, the only relevant 
terms that are lacking from the ab initio description are 
the polarization terms. In a polarizablc point-ion approx- 
imation, the polarization contribution to the deformation 
and interaction energies is 
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where is the polarizabihty of the ion R, and pin the 
dipole moment induced on ion R. The new terms added 
to the interaction energy are the monopole-dipole and 
dipole-dipole interaction energy terms. The term added 
to the deformation energy represents the energy cost of 
deforming the charge density of the ion to create the 
dipole moment. The point-ion approximation provides 
just the asymptotic part of the polarization interaction 
energy, that is, it is exact only for large ionic separa- 
tions. As soon as the ions begin to overlap, there is an 
important-short-range contribution to the induced dipole 
moments, circ3 which is of opposite sign to the asymptotic 
limit for anions and may in some specific cases reverse 
the sign of the asymptotic value. These effects can be 
easily acomodated in the formalism by substituting the 
asymptotic value of the induced dipole moments by the 
following expression: 
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where sr stands for "short-range" and ^^o*'''-^ is the a 
component of the dipole moment vector induced on ion 
R. The physics behind the short-range polarization cor- 
rection has been explained in Ref. p]j4^ , and is associated 
to the finite extents of the electron densities of the an- 
ions and cations. Then, / is a short-range function that 
switchs on as the cation-aniop-overlap becomes apprecia- 
ble. Madden and coworkeraij have emp loved the Tang 
and Toennies dispersion damping functioned as a suitable 
form for /: 
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This is a smoothed step function passing from zero for 
large r to — c for r=0. The range of r values at which 
/ becomes significantly different from zero is primarily 
determined by the range parameter b. 

We have included the polarization terms in the en- 
ergy calculation with this parameterised method, that 



calculates the induced dipole moments from eq. (5) and 
the correction to the deformation and interaction ener- 
gies from eqs. (3) and (4), respectively. The "enlarged" 
aiPI-|-polarization model thus obtained accounts for all 
the relevant physical interactions. The relaxation of the 
assumption of spherical symmetry being computation- 
ally expensive, the price to be paid is the inclusion in the 
model of a set of parameters, namely, the polarizabili- 
ties aji and the range parameter b. Appropriate values 
for the other two constants c and kmax caa be taken 
equal to the bulk values (c=-3 and kmax—4:).LB Given the 
meaning of the range parameter 5, inversely related to 
half the interionic distance between first neighbors, one 
might expect different values of b for clusters as com- 
pared to bulk materials if the interionic distances are 
substantially different. As a matter of fact, the evolu- 
tion of those distances with cluster size-ia_not too com- 
plicated in the case of ionic clusters B3'E£rM Specifically, 
the average interionic distance d initially increases quite 
abruptly with the number of molecules n, and then slowly 
approaches the bulk limit. As a consequence of this bet. 
haviour, we will see that the bulk value (6=0.75 a.u.)c£l 
is appropriate for all (CaO)„Ca^+ clusters with n > 4. 
Different values of b are needed just for n <4 to avoid 
overpolarization problems.cB Regarding the polarizabili- 
ties, oxide anions have the interesting property of show- 
ing strongly coordination-dependent values. In fact, the 
O'^" anion does not exist as a free ion, which is equiv- 
alent to an infinite polarizabihty; in the solid phase it 
is stabilized by the crystal environment £pd has a finite 
material-dependent polarizabihty. Wilsor£3 has interpo- 
lated between those two limits and gives values for the 
coordination-dependent values of a{0^~) in MgO. We 
have assumed that the ratio of the bulk oxide polarizabil- 
ities for MgOj^d CaO (^''""^(O^- ; CaO)/a''""=(02- : 
M(70)=1.469)cJ is independent of the oxide coordina- 
tion, and have deduced the a values for CaO from those 
of MgO. This procedure is justified because the cation 
size does not change appreciably with coordination num- 
ber. For the caicium cation we take the bulk polarizabil- 
ity (3.193 a.u.)El 

We close this section with a consideration of several 
criticisms that could be raised against (and of the ad- 
vantages of) the employed methodology. We have cho- 
sen a mixed ah mitio/semiempirical energy model in or- 
der to obtain a good compromise between computational 
efficiency and accuracy. All the relevant energy terms 
excluding polarization are described with an ab initio 
methodology. To. include polarization, we have used an 
accurate model,c3 where the parameters haye been fitted 
by a comparison to ab initio calculations £3 Special care 
has been devoted to the separation of all the indepen- 
dent physical factors that influence a given quantity, thus 
avoiding a mixing of different effects in a single parameter 
and enhancing the transferability of the model. The good 
paramctcrisation is reflected in the fact that parameters 
can be transfered between closely related systems (like, 
for example, different metal oxides) by simple scaling ar- 
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guments involving ionic radii.ca Thus, we think that the 
rehabihty of our calculations is reduced just a little com- 
pared to full ab initio methodologies. To support this 
expectation, we made a comparison with DMOL calcu- 
lations performed as. neutral (CaO)„ clusters by Malli- 
avin and Coudray.cJ All the interionic distances were 
in agreement to their calculations up to differences of 4 
%. The energetic ordering of the isomers, as well as the 
specific energy differences, are reproduced with a maxi- 
mum error of 5 %. We believe that this is a very rea- 
sonable agreement, even more if we realize that we are 
neglecting dispersion interactions, and polarization in- 
teractions beyond the dipolar terms. The solid MgO is 
excellently described with the aiPI model (at least in its 
static properties) .ca The model is then expected to trans- 
fer properly between both limits. The larger computa- 
tional simplicity has been exploited to study large cluster 
sizes (up to 59 ions) with full relaxations of the geome- 
tries. Moreover, for each cluster size, a large number 
of isomers (between 10 and 15) have been investigated. 
The generation of the initial cluster geometries was ac- 
complished by using a pait. potential, as we explained in 
our previous publication.La The optimization of the ge- 
ometries haspbeen performed by using a downhill simplex 
algorithm.ElEl 

III. RESULTS AND DISCUSSION 

A. Structural Trends in (CaO)„Ca'^+ Cluster Ions 

In Fig.l we present the optimized aiPI+polarization 
structures of the ground state (GS) and lowest lying iso- 
mers or (CaO)„Ca^+ (n=4-29) cluster ions. Below each 
isomer we show the energy difference (in eV) with re- 
spect to the ground state. For n<4, the clusters are not 
detected in the experiments, probably because they un- 
dergo a Coulomb explosion driven by the excess charge 
and the small cluster size, but we are not interested here 
in this aspect of the experiments. From n=4 to n=10, 
there is a predominance of mx2x2 fragments (m=2~5), 
that is, the (CaO)2 subunit appears as the basic building 
block. The total number of ions in these nonstoichiomet- 
ric clusters is an odd number, and thus those structures 
are never perfectly compact. There is either an extra 
cation added to or a missing anion removed from the 
perfect structure. Less compact structures as for example 
planar fragments are not energetically competitive. The 
structures in this size range tend to be elongated as a di- 
rect consequence of the excess cluster charge. When n=9, 
a 3x3x2-1-1 fragment is more stable than that based on 
the (CaO)2 building block, and n=10 is the largest clus- 
ter size for which a fragment of this kind is the ground 
state. Another thing to be pointed out is that in this size 
range, the extra cation present in the mx2x2-|-l struc- 
tures induces a larger cluster distortion than the missing 
anion in the mx 2x2-1 structures. We will see that this 



feature has important implications in the stability of the 
clusters. 

From n=ll to n=:15, the dominant fragments are based 
on mx3x2 units. For n=16 and 17, the most stable iso- 
mers are mx4x2 fragments. None of these structures has 
still developed an anion with full bulk coordination. In 
particular, the 3x3x3 isomer for n=13, which is particu- 
larly stable in. the case of nonstoichiometric alkali halide 
cluster ions,Eil does not even appear in Fig. 1. The large 
coordination-dependent values of the polarizabilities of 
the oxide anions favors the formation of surface sites, 
and gives rise to somewhat less compact ground state 
structures, for which the increase in dipolar energy com- 
pensate for the decrease in Madelung energy. The 3x3x3 
(CaO)i3Ca^+ is specially unfavored by the dipolar energy 
terms because it has a central oxide anion with bulk co- 
ordination (so with a comparatively low polarizability) , 
and another 12 anions with coordination 4. On the con- 
trary, the largest coordination in the GS structure is five, 
and some three-coordinated anions (in corner positions) 
also appear, inducing a large dipolar energy stabilization. 
For n=18 and 19 there is a glimpse of a transition to more 
compact cluster structures. The important feature of the 
GS structures of these two cluster sizes, compared to the 
3x3x3 for n=13, is that now there are oxide anions in 
corner positions. These make a large contribution to the 
polarization energy term, that added to the increased 
Madelung energy of a compact fragment, gives a total 
GS energy more negative than that ofmx3x2ormx4x2 
structures. Nevertheless, the energy differences between 
isomers are small, and from n=20 on, ground state iso- 
mers without bulk anions are again obtained (n=24 and 
27 are the only relevant exceptions, because the ground 
states of n=26 and n=29 can be considered degenerate 
within the accuracy of our theoretical model). 

A general feature of (CaO)„Ca^"'" cluster ions with n>8 
is that axbxc-f 1 fragments are specially stable compared 
to other isomers whenever they can be formed. In Ta- 
ble I we show all the fragments of that kind relevant to 
the cluster size range considered in this study. Each se- 
ries has a typical periodicity that could in principle be 
refiected in different portions of the mass spectra, given 
the high stability of these fragments. Some sizes can be 
accomodated in several families, that is, the classification 
is highly redundant, but useful anyway to our purposes. 
If for a given cluster size, a cluster with that formula can 
be formed, it is always the ground state structure. If it 
is possible to build up two different isomers with that 
formula (n=12, 18, 24), the more compact structure is 
energetically favored. This rule works as long as we do 
not consider structures that are not energetically compet- 
itive anymore (the isomer based on the (CaO)2 building 
block of (CaO)i4Ca^+ is an example), and can be helpful 
in guessing specially stable structures for clusters larger 
than those studied here. For nearly all those cluster sizes 
with no competitive axbxc+1 structure, axbxc-1 frag- 
ments are obtained as the ground state or specially stable 
isomers (examples are found for n=5, 7, 11, 19, 23 and 
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29). The special stability of axbxc+1 structures is some- 
times reflected in high stabilities for the corresponding 
axbxc-|-3 structures, comparable indeed to the stabili- 
ties of axbxc-1 fragments; this occurs for n=13, 17, 19 
and 26. With the only exceptions of n=13,14,22,26 and 
29, all (CaO)„Ca^+ GS structures are explained in terms 
of those three kinds of fragments. 

Comparing to the results, of our previous paper on 
(MgO)„Mg^+ cluster ions,L2l we can see that from n=4 
to n=20 the GS structures are basically the same in both 
systems (the only exceptions are ii—7 and n=13). Inter- 
esting structural differences between both materials ap- 
pear in the size range n>20. Specifically, the transition 
to bulklike structures, containing inner anions with bulk 
coordination, is slower in the case of (CaO)„Ca^+. An 
analysis of the several energy components shows that the 
net effect of polarization is more important in calcium 
oxide. Although the polarizability of Ca^^ is larger than 
that of Mg^+ and the coordination-dependent values of 
a(0^~) are larger in CaO than in MgO, this is not a triv- 
ial conclusion, because the interatomic distances are also 
larger in CaO, and so the electric fields acting on each ion 
are correspondingly smaller. If we consider that a highly 
ionic material is that one for which the Madelung en- 
ergy term is almost completely dominant in determining 
structural and several other properties, we would con- 
clude that the ionic character of (CaO)„Ca^"'" is smaller 
that that found for (MgO)„Mg^+. Indeed, being the po- 
larization contribution more important, the structures of 
calcium oxide clusters have a larger directionality degree, 
a feature that is usually associated to covalency (opposite 
to the natural tendency of purely ionic systems to form 
isotropic structures). However, we think that the term 
"covalency" should be employed just in those situations 
where charge transfer between different atomic centerSjia 
important. As Madden and coworkers have discussed,C^ 
polarization terms in ionic systems are responsible for a 
lot of properties traditionally attributed to "covalency" . 
One point that deserves further investigation, however, 
is whether the directional properties induced by polar- 
ization effects (and reflected in a lower average coordi- 
nation) can be responsible for a larger charge transfer 
between different centers. This would be reasonable be- 
cause the saturation of the bonds is less complete. 

B. Relative stabilities and connection to 
experimental mass spectra 

In the experimental mass spectra, the populations 
observed for some cluster sizes are enhanced over those 
of the neighboring sizes. These "magic numbers" are a 
consequence of the evaporation events-that occur in the 
cluster beam, mostly after ionization^ A magic cluster 
of size n has an evaporation energy that is large compared 
to that of the neighboring sizes (n-1) and (n-|-l). Thus, 
on the average, clusters of size n undergo a smaller num- 
ber of evaporation events and this leads to the maxima 



in the mass spectra. As our main concern in this section 
is to compare with the experimental mass spectra, we 
calculate the evaporation energy as a function of cluster 
size. To do this, we assume that the dominant evap- 
oration channel is the loss of a neutral (CaO) molecule, 
something supported by the experiments of Ziemann and 
Castleman.lia In the size range n<ll, some other chanels 
seem to be opened in the experiments,t3 and indeed for 
n<4 Coulomb explosion is dominant, that is the reason 
why we do not consider clusters with n<4. With that as- 
sumption, the evaporation energy of (CaO)„Ca^+ reads 

Eevapin) = Eel uster 

[{CaO)n-iCa^+]+E{CaO) 

-E,iu,teAiCaO)nCa^+]. (9) 

Maxima in the evaporation energy curve do not al- 
ways coincide with maxima in the experimental mass 
spectraEa There are two main processes that contribute 
to enhance the cluster population for size n: a)A small 
evaporation energy for size (n+1); b)A large evaporation 
energy for size n. Thus, a most convenient quantity to 
compare with experiment is the second energy difference 

A2(n) = Eavapin + 1) - Eevapin). (10) 

A negative value of A2 (n) indicates that the n-population 
increases by evaporations from the (n-l-l)-clusters more 
rapidly than it decays by evaporation to the (n-1)- 
clusters. Specifically, the specially stable cluster sizes 
will be reflected as minima in the A2(n) curve. 

Now, the evaporation energy Eevap{n) of eq. (9) can 
be calculated in two different ways. In the first one, 
energy differences are always taken between the ground 
state structures of sizes n and (n-1). This procedure, 
which we call (by obvious reasons) adiabatic evapora- 
tion, reflects the stability of the clusters in the limit of 
small energy barriers between isomers or alternatively of 
large experimental times of flight. The stabilities calcu- 
lated in this way are shown in the upper part of figure 
2. Magic numbers are found for n=5,8, 12, 15, 18,20,24,27. 
The only axbxc-1 structure that shows a special stabil- 
ity is that of n=5. The rest of magic clusters belong to 
the axbxc-|-l family of structures. If n is a magic size, 
and both (n-l-l) and (n-1) GS structures do not belong 
to the axbxc+1 family, a deep minima is found in the 
A2{n) curve (this happens for n=12 and 18). For the 
rest of magic sizes, the (n-|-l) GS structure has also the 
formula axbxc+1, and has a correspondingly high sta- 
bility reflected in a negative value of A2(n + 1). In these 
cases the stability of size n is just slightly enhanced over 
that of size (n+1). One can appreciate the increasing rel- 
evance of the Madelung term in determining the cluster 
stabilities: when n<20, the most stable axbxc+1 struc- 
tures are the less compact ones (n=8 and 15 more stable 
than n=9 and 16, respectively); if n>20, that trend is 
reversed (n=20, 24 and 27 more stable than n= 21, 25, 
and 28, respectively). The special relevance of axbxc+1 
structures in explaining the cluster stabilities does not 
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conform to the initial experimentaLexpectations of high 
stabihties for axbxc-1 structuresJij Analysing the en- 
ergy components, we find that the polarization contri- 
bution stabilizes the axbxc+1 structm^e more than the 
corresponding axbxc-1 structure for all values of a,b,c. 
For the smallest cluster sizes, however, the extra cation 
present in axbxc+1 structures induces a large cluster 
distortion compared to that induced by the missing anion 
in axbxc-1 structures, and the Madelung contribution 
favors these last structures in a larger amount, making 
them more stable for some sizes. 

The second kind of calculation of Ee„ap(n) proceeds 
as follows: we consider the optimized GS structure 
of (CaO)„Ca^+ and identify the CaO molecule that 
contributes the least to the cluster binding energy. 
Then we remove that molecule and relax the resulting 
(CaO)„_iCa^+ fragment to the nearest local minimum. 
This process can be termed locally adiabatic because 
both fragments are allowed to relax to the local min- 
imum energy configuration after the evaporation. For 
some cluster sizes, the fragment of size (n-1) left when 
a CaO molecule is removed from (CaO)„Ca^'*' does not 
lie on the catchment basin of the (CaO)„_iCa'^+ GS iso- 
mer, so that the locally adiabatic evaporation energies 
are larger than the energy differences between adjacent 
ground states minus E(CaO) in those cases. The locally 
adiabatic evaporation energies are plotted as a function 
of n in the lower part of Fig. 2. These will reflect the clus- 
ter stabilities in the limit of large energy barriers between 
isomers or of short experimental times of flight. Magic 
numbers are obtained for n=5, 7, 9, 11, 13, 16, 19, 22, 25 and 
27, in complete agreement with the experiments of Zie- 
mann and Castleman.Ej 

The main message to be extracted from these consider- 
ations is that the magic numbers obtained in the experi- 
ments might be dominated by the effects of kinetic traps 
occuring in the course of the evaporation process. Since 
our calculations are static, we can not rigorously assert 
that this is the only possible explanation, but a plausibil- 
ity argument based on a comparison to the closely related 
and more thoroughly studied case of alkali halides sup- 
ports our expectations. The mpibility experiments per- 
formed by the group of JarroldaEl show that the relax- 
ation dynamics to the ground state structure for sodium 
chloride clusters involves drift times of almost one sec- 
ond. The importance of kinetic traps in explaining these 
interesting resultsL.is shown in the theoretical works of 
Doye and Wales O Specifically, these authors show that 
the potential energy landscape of alkali halide clusters, 
calculated by using a phenomenological pair potential to 
describe the interactions, is structured in several funnels, 
separated from each other by high free-energy barriers. 
When a cluster evaporates a molecule, it cools in the 
process, so trapping kinetic effects are expected when- 
ever parent and product GS structures belong to different 
funnels. Given the close similarities between halide and 
oxide systems, one expects similar effects in the evap- 
oration kinetics of alkaline-earth oxides to be relevant. 



The main structural differences are due to the effects of 
polarization, and these could also affect the mechanisms 
of structural transitions. In the case of alkali halides, 
Doye and Wales find that a highly cooperative process 
is energetically less impeded by energy barriers than se- 
quential ionic diffusion,Ej with interesting implications 
for the mechanical properties of these clusters. Perhaps 
the same is true for the clusters studied here, but one 
has to keep in mind that polarization tends to lower the 
barriers against diffusion,E3 and those effects are more 
important for oxides. We think that further calculations 
of this kind for oxide clusters would be very interesting. 
Mobility experiments on (MgO)„Mg^+ or (CaO)„Ca^+ 
could conclusively confirm the structural trends found in 
the present work. 

C. Neutral Stoichiometric (MgO)„ and (CaO)„ 
clusters 

The experiments performed by Saundersi0 show that 
both (MgO)+ and (CaO)+ stoichiometric cluster ions 
with a number of molecules n=6,9,12 and 15 are expe- 
cially abundant in the mass spectra. However, when 
these clusters are allowed to collide with inert gas 
ions, the fragmentation channels are different: (MgO)3 
fragments are predominantly observed in one case and 
(CaO)2 fragments in the other. These results suggest 
that the basic cluster building blocks are different for the 
two materials, but not so different as to lead to different 
magic numbers. 

We foLLud a similar scenario in the case of alkali halide 
clusters.EHl Specifically, a universal set of magic numbers 
n=4,6,9,12,... was found for the whole family of (AX)„ 
clusters, with A=Li,Na,K,Rb and X=F,Cl,Br,I. However, 
the cluster structures were not found to be the same 
for all the different materials. When the cation size is 
much smaller thMi the anion size (all lithium halides and 
sodium iodide) J23 ground state structures based on the 
stacking of hexagonal (AX)3 rings are obtained. For 
the rest of materials, the ground state structures are 
mostly obtained by stacking of rectangular (or double- 
chain) (AX)3 planar fragments. This is just a packing ef- 
fect: when the ratio of cation to anion size is very small, 
anion-anion overlap repulsive interactions are large, forc- 
ing an opening of the (AX) 3 rectangular fragments into 
hexagons. The magic numbers are the same for both 
structural families because it is for those cluster sizes that 
specially compact structures can be formed,-. When we 
studied (AX)„A"'" alkali halide cluster ionsj^l we found 
that the structures were much more similar irrespective 
of packing considerations. The ring structures are not 
competitive in this case because it is not possible to build 
up a perfect hexagonal fragment with an odd number of 
ions. On the contrary, perfect cubic structures can be 
formed (as for example the 3x3x3 structure for n=13). 

From our study on doubly-charged clusters, we have 
not found important structural differences between 
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(CaO)„Ca2+ and (MgO)„Mg2+,ii at least in the small 
size regime. We have performed aiPI+polarization cal- 
culations on the structures of (MgO)„ and (CaO)„ with 
n=3,6,9,12,15 and 18. Specifically, we have considered 
just those structures based on staking of hexagonal and 
rectangular (A0)3 units, and those based on stacking of 
(A0)2 units, with A=Mg or Ca. We find that the ground 
state structures of (MgO)„ clusters are based on (MgO)3 
units, while those of (CaO)„ clusters have a rectangular 
(CaO)3 building block, being this the same packing ef- 
fect found in the case of alkali halides. The structure of 
(CaO)6 could be alternatively viewed as the stacking of 
three (CaO)2 units, but for n=9,12,15 and 18, the tubu- 
lar shapes obtained by stacking (CaO)2 units are not 
competitive anymore. Were all the ground state struc- 
tures of (CaO)„ clusters based on the (CaO)2 building 
block, we would expect a periodicity of 2 in the magic 
numbers observed in the mass spectra. Saunders shows 
the collision induced fragmentation spectra of (CaO)„, 
with n=4,6,8,Ll which are certainly based on stacking of 
(CaO)2 units, but does not show those for (CaO)9 or 
(CaO)i2, for example. Our main conclusion is that the 
special stability of (CaO)„ clusters is also explained in 
terms of (CaO)3 units, but with rectangular instead of 
hexagonal shape. This explains the same periodicities 
observed in the magic numbers of both materials. 

IV. SUMMARY 

The ab initio perturbed ion model, supplemented with 
a parameterised treatment of dipolar terms, has been 
employed in order to study the structural and energetic 
properties of (CaO)„Ca^"'" (n=l-29) cluster ions. Po- 
larization effects favor the formation of surface sites, and 
reduce the stability of highly compact structures contain- 
ing anions with bulk coordination. Thus, despite many 
similarities in the experimental mass spectra, the struc- 
tures of alkaline-earth oxide and alkali halide cluster ions 
are shown to be different. Most of the lowest energy 
structures have the formula axbxc-l-1. The structures 
of (CaO)„Ca^'*' and (MgO)„Mg^+ cluster ions are very 
similar for n<20, irrespective of differences in cationic 
size and polarization. It is just for n>20 that structural 
differences emerge, showing a slower convergence to bulk 
properties for CaO compared to MgO. The analysis of 
the stabilities suggests that the experimental mass spec- 
tra could be dominated by the effects of kinetic traps. 
Specifically, if we consider locally adiabatic evaporation 
events, complete agreement is found with the experi- 
mental stabilities. The neutral stoichiometric (MgO)„ 
and (CaO)„ clusters (n=3,6,9,12,15,18) show structural 
differences similar to those observed in neutral stoichio- 
metric alkali halide clusters: the basic building block is 
an (MgO)3 hexagonal fragment in the case of MgO and 
a (CaO)3 rectangular (or double-chain) fragment in the 
case of CaO. This is just a packing effect due to the 



larger overlap repulsion between anions when the cation 
size is very small. While the structures of (CaO)„ clus- 
ters, with n=4,6,8 are certainly based on (CaO)2 units, 
as suggested by collision-induced fragmentation experi- 
ments, the specially stable (CaO)„ clusters are based on 
a (CaO)3 unit. This explains the same periodicity of 
3 observed in the experimental magic numbers of both 
(MgO);^ and (CaO)J clusters. 

Captions of Figures and Tables. 

Figure 1. Lowest-energy structure and low- lying iso- 
mers of (CaO)„Ca^"'" cluster ions. Dark balls are Ca^"*" 
cations and light balls are 0^~ anions. The energy dif- 
ference (in eV) with respect to the most stable structure 
is given below the corresponding isomers. 

Figure 2. Adiabatic (a) and locally adiabatic (b) 
evaporation energies required to remove a neutral CaO 
molecule from (CaO)„Ca^+ cluster ions as a function of 
n. The local minima in the evaporation energy curve are 
shown explicitely. 

Table I Possible different axbxc-|-l structures, with 
their inherent periodicities. Those cluster sizes n that 
are actually observed as ground state structures of 
(CaO)„Ca^+ clusters are written in boldface. 



Structure 


Periodicity 


Cluster size n 


mx2x2-|-l 


2 


8,10,12,... 


mx3x2-|-l 


3 


9,12,15,18,21,24,27... 


mx4x2-|-l 


4 


8,12,16,20,24,28,... 


mx5x2-|-l 


5 


10,15,20,25,... 


mx6x2+l 


6 


12,18,24,... 


mx3x3+l 


9 


9,18,27,... 


mx4x3+l 


6 


12,18,24,... 
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